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We report a quark spin calculation from the anomalous Ward identity with overlap fermions on 
2+1 flavor dynamical fermion configurations with light sea quark masses. Such a formulation de¬ 
composes the divergence of the flavor-singlet axial-vector current into a quark pseudoscalar term 
and a triangle anomaly term, flavor by flavor. A large negative contribution from the anomaly 
term is observed and it is canceled within errors by the contribution from the pseudoscalar term 
in the disconnected insertion in the heavy quark region. On the other hand, net negative contribu¬ 
tions are obtained for the light and strange quarks in the disconnected insertion, since their quark 
pseudoscalar terms are smaller than that of the heavy quark. Our results are obtained from the 
2+1 flavor domain wall fermion configurations on the 24 3 x 64 lattice with a -1 = 1.78(5) GeV 
and the light sea quark at m n = 330 MeV. We use the overlap fermion for the valence and the 
quark loop so that the renormalization constants Z m and Zp cancel in the pseudoscalar operator 
2 mP. In addition, the overlap Dirac operator is used to calculate the local topological charge in 
the anomaly so that there is no renormalization for the anomaly term either. In this study, we find 
the total quark spin to be small mainlyly due to the large negative anomaly term which could be 
the source for the ‘proton spin crisis’. 
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1. Introduction 


Apportioning the spin of the nucleon among its constituents of quarks and gluons is one of the 
most challenging issues in QCD both experimentally and theoretically. It is shown by X. Ji [9] that 
there is a gauge-invariant separation of the proton spin operator into the quark spin, quark orbital 
angular momentum, and glue angular momentum operators 


J QCD — Jq~\~ Jg — 2 "3" Lq T Jgi 


( 1 - 1 ) 


where the quark and glue angular momentum operators are defined from the symmetric energy- 
momentum tensor 


Ji 


q,g 




(1.2) 


with the explicit expression 


Jq — ^ + L q — f d X 




(1.3) 


for the quark angular momentum which is the sum of quark spin and orbital angular momentum, 
and each of which is gauge invariant. 

Since the contribution from the quark spin is found to be small (~25% of the total proton 
spin) from the global analysis of deep inelastic scattering data [1], it is expected that the remain¬ 
der should come from glue spin and the orbital angular momenta of quarks and glue. The quark 
spin contribution from u, d and s has been studied on the lattice [2, 3] since 1995 with quenched 
approximation or with heavy dynamical fermions [4], Recently, it has been carried out with light 
dynamical fermions [5, 6, 7, 8] for the strange quark. The strange quark spin (As) they found, is in 
the range from —0.02 to —0.03 which is several times smaller than that from a global fit of DIS and 
semi-inclusive DIS (SIDIS) which gives As ss —0.11 [1]. Such a discrepancy between the global 
fit of experiments and the lattice calculation of the quark spin from the axial-vector current has 
raised a concern that the renormalization constant for the flavor-singlet axial-vector current could 
be substantially different from that of the isovector axial-vector current [10, 11] since the latter is 
commonly use for the lattice calculations of the flavor-singlet axial-vector current for the quark 
spin. To alleviate this concern, we use the anomalous Ward identity (AWI) to calculate the quark 
spin. In this proceeding, we report the calculation of both the connected insertion (Cl) and discon¬ 
nected insertion (DI) contributions to quark spin from u, d, s and c, which are non-perturbatively 
renormalized when the anomalous Ward identity is used to carry out the calculation with the over¬ 
lap fermion for the quarks and the overlap Dirac operator for the topological charge. 

The anomalous Ward identity includes the triangle anomaly in the divergence of the flavor- 
singlet axial-vector current 

Nf 

= 2 £ mfq f iY$q f + 2iNfq, (1.4) 

/=i 

where q is the local topological charge operator and is equal to j^Jr c G^ v G^ v in the continuum. 
We put this identity between the nucleon states and calculate the matrix element on the right-hand 
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side with a momentum transfer q 

jM 

(p's lAjj ps)sn = lim Qr\(p',s\2 Y m f qfiy 5 q f + 2iN f q\p,s). (1.5) 

9->o q ■ s “| 

Lattice QCD theory has finally accommodated chiral symmetry, the lack of which has hampered the 
development of chiral fermions on the lattice for many years. It is shown that when the lattice mass¬ 
less Dirac operator satisfies the Gingparg-Wilson relation y$D + Dy$ = aDy$D, for which the over¬ 
lap fermion with negative mass parameter being the explicit example [14], the modified chiral trans¬ 
formation leaves the action invariant and gives rise to a chiral Jacobian factor J = e~ 2iaTr/si l ~2 aD > 
from the fermion determinant [12]. The index theorem [13] shows that this Jacobian factor carries 
the correct chiral anomaly. It is shown further that the local version of the overlap Dirac operator 
gives the topological charge density operator in the continuum [15], i.e. 

7>y 5 (l - * aD ov (x,x )) = j^tr c G^ v G^ v (x) + @(a 2 ) (1.6) 

Therefore, Eq. (1.4) is exact on the lattice for the overlap fermion which gives the correct anoma¬ 
lous Ward identity at the continuum limit. Instead of calculating the matrix element of the axial- 
vector current derived from the Noether procedure [16, 13], we shall calculate it from the r.h.s. of 
the AWI in Eq. (1.4) through the form factors defined in Eq. (1.5). 

In the lattice calculation with the overlap fermion, we note that the renormalization constant of 
the pseudoscalar density cancels that of the renormalization of the quark mass, i.e. Z m Zp = 1 for the 
chiral fermion. Also, the topological charge density, when calculated with the overlap operator as 
in the l.h.s of Eq. (1.6) is renormalized - its integral over the lattice volume is an integer satisfying 
the Atiya-Singer theorem. Thus, when the matrix elements on the right-hand side of Eq. (1.5) are 
calculated with the overlap fermion and its Dirac operator, the flavor-singlet axial-vector current is 
automatically renormalized on the lattice non-perturbatively ci la anomalous Ward identity (AWI). 

2. Numerical details 

Besides the fact that AWI admits non-perturbative renormalization on the lattice, the pseu¬ 
doscalar density and the topological density represent the low-frequency and high-frequency parts 
of the axial-vector quark loop respectively. On the 2+1 flavor domain-wall fermion configurations 
on the 24 3 x 64 lattice (a -1 = 1.78(5) GeV [17]) with the light sea quark mass corresponding to 
a pion mass at 330 MeV. It is learned that on such a lattice, a mere 20 pairs of the overlap low 
eigenmodes would saturate more than 90% of the pseudoscalar loop with light quarks in config¬ 
urations with zero modes [18]. On the other hand, it is well-known that the contribution to the 
triangle anomaly comes mainly from the cut-off part of the regulator. Therefore, the topological 
charge density represents the high-frequency contribution of the axial-vector loop, albeit in a local 
form (the overlap operator is exponentially local). Since the quark loop of the pseudoscalar density 
is totally dominated by the low modes, we expect that our low-mode averaging (LMA) approach 
should be adequate for this term. To the extent that the signal for the topological term is good, one 
should be able to calculate the flavor-singlet g a with the AWI . Both the overlap fermion for the 
quark loop and the overlap operator for the topological charge density are crucial in this approach. 
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Figure 1: The summed matrix elements for the light quark loop (corresponds to m K = 330 MeV) with the 
pseudoscalr current, as a function of the separation of the source and sink. We confirm that the low mode 
part dominates in the pseudoscalar case. 

We shall first check to see how the low mode averaging works for quark loops of the pseu¬ 
doscalar operator. We plot the ratio of the 3-pt to 2-pt correlation functions in the DI for the sum 
method where the insertion is summed between the nucleon source and sink times so that the ma¬ 
trix element is the slope of the ratio at large source-sink time separation. This is done for the case 
where both the quarks in the loop and in the nucleon coincide with that of the light sea mass which 
corresponds to m n = 330 MeV. 

The reason for the good signal is that, besides adopting LMA in the DI, that low-mode substi¬ 
tution (LMS) technique with 8 smeared noise-grid source on a given time slice [18] has been used 
to improved the nucleon propagator in the construction of the three-point function, and all the time 
slices has been looped over to improve the statistics. 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 
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Strange (17^=330 MeV) 



(a) (b) 

Figure 2: (a) The charm pseudoscalar and topological density contributions to the proton spin as a function 
of Q 2 . (b) The same as in (a) for the strange. 

We first show the results for the charm quark which contribute only in the DI. The pseudoscalar 
density term and the topological charge density term are plotted in Fig. 2a as a function of Q 2 . We 
see that the pseudoscalar is large due to the large charm mass and positive, while the topological 
charge term is large and negative. When they are added together (black triangles in the figure), it is 
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consistent with zero for the whole range of Q 2 . When extrapolated to Q 2 = 0, the charm gives zero 
contribution to the proton spin within error due to the cancellation between the pseudoscalar term 
and the topological term. It is shown [19] that the leading term in the heavy quark expansion of the 
quark loop of the pseudoscalar term, i.e. 2 mP is the topological charge tr c G^ v G^ v , but with a 
negative sign. Thus, one expects that there is no contribution to the quark spin from heavy quarks 
to leading order. It appears that the charm quark is heavy enough so that the G(\/m 2 ) correction is 
small. We take this as a cross check of the validity of our numerical estimate of the DI calculation 
of the quark loop as well as the anomaly contribution. 

The contributions from the strange are also calculated and shown in Fig. 2b. The 2 mP con¬ 
tribution is slightly smaller than that of 2 q and results in a net small negative value for the sum of 
2 mp and 2 q at finite Q 2 . After a dipole fit, we obtain A.s = —0.26(5) at m n = 330 MeV. 

m q =500 MeV light iso-vector (1^=330 MeV) 




(a) (b) 

Figure 3: (a) The quark spin combinations of the proton-like baryon with m q ~ 500 MeV, in which the DI 
contribution is canceled by the topological change one. (b) The same as in (a) for light quarks at the unitary 
point. 

Since this A.s is quite a bit smaller than the experimental value, we explore the possible finite 
volume effect and the fact that the induced pseudoscalar form factor /?a(^ 2 ) has been neglected 
in the Q 2 extrapolation which does not contribute at the Q 2 = 0 limit as in Eq. (1.5), but has a 
contribution at finite Q 2 [20]. We shall check this in the connected insertion (Cl) calculation. As 
can be seen in Fig. 3 for m q ~ 500 MeV, both A u and Ad in Cl calculated from the axial-vector 
current and renormalized with Z A from the isovector Ward identity are well reproduced through the 
Q 2 extrapolation of 2 mP with a dipole form. Whereas, in the case of light quarks at the unitary 
point, = 1.13(2) from the axial-vector current is 1.8(1) times larger than 0.62(4) from the 
dipole extrapolation of 2mP. This is most likely due to our ignoring of the induced pseudoscalar 
form factor li A (tj 2 ) as well as the finite volume effect at small Q 2 which is well known to plague 
the extrapolation of nucleon magnetic form factor. 

At the unitary point, when the valence u/d mass matches with that of the light sea, Fig. 2 
shows the combined pseudoscalar contributions 2 mP u / d from both the combined Cl and DI with a 
dipole extrapolation. Also plotted are the overall quark spin Au/Ad by including the topological 
charge contribution. In this case, we obtain Au + Ad = 0.19(3) and Au — Ad = 0.62(4) at Q 2 = 0 
from a dipole extrapolation in Q 2 . As we discussed above, the fact that from the axial current 
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light (01^=330 MeV) 



Figure 4: The combined pseudoscalar contribution from both the connected insertion (Cl) and DI (2 mP u / d 
in the plot), along with the overall quark spin from both pseudoscalar and topological charge (g A ,d)- The plot 
corresponds to the unitary point with m n = 330 MeV. 

is 1.8(1) times larger than that of An — Ad through the Ward identity approach is most likely due 
to the neglect of the induced pseudoscalar form factor /?a(^ 2 ) and the finite volume effect in the 
Q 1 extrapolation. We apply this 1.8(1) factor as an estimate to correct the present AWI approach 
and obtain Au + Ad = 0.35(6), As = —0.05(1). Thus the total estimated spin A£ = 0.30(6) at the 
unitary point is consistent with the present experimental results which are between 0.2 and 0.3. We 
expect that, at lighter quark masses, AT will be smaller. 

The above results are from the 24 3 x 64 lattice with 200 configurations. The nucleon propa¬ 
gator in the DI has been calculated with the smeared-grid noise source with time dilution which 
covers all time slices in order to have reasonable statistics for the DI. 

3. Summary 

We have carried out a quark spin calculation from the anomalous Ward identity with overlap 
fermions on 2 + 1 flavor dynamical fermion configurations with light sea quark masses. A large 
negative contribution from the anomaly term is observed and it is canceled within errors by the 
contribution from the pseudoscalar term in the heavy quark region. On the other hand, net negative 
contributions are obtained for the light and strange quarks in the disconnected insertion, since their 
quark pseudoscalar terms are smaller than that of the heavy quark. Since the overlap fermion is used 
for the pseudoscalar term 2 mP and the overlap Dirac operator is used for the local topological term, 
the result from AWI is non-perturbatively renormalized. Since in the Cl, the dipole extrapolation 
of 2mP obtained a An — Ad which is smaller than that obtained from axial vector current with 
renormalization Z A from the Ward identity, one is concerned that there is substantial finite volume 
effect and also the fact that the induced pseudoscalar form factor li A (q 2 ) is neglected in the Q 2 
extrapolation. These effects need to be studied in the upcoming calculation on the 48 4 X 96 lattice 
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which has twice the spatial size (i.e. 5.5 fm) of the present one and the pion mass is at the physical 
value. 

As far as the smallness of the quark spin in the proton is concerned, present results from only 
one lattice notwithstanding, has presented hints to suggest that the culprit of the ‘proton spin crisis’ 
is the 1/(1) triangle anomaly. 
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